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We study the nonlinear dynamics of an ensemble of cold trapped atoms with a hyperfine transition 
magnetically coupled to a resonant microwave cavity mode. Despite the minute single atom coupling 
one obtains strong coupling between collective hyperfine qubits and microwave photons enabling 
coherent transfer of an excitation between the long lived atomic qubit state and the mode. Evidence 
of strong coupling can be obtained from the cavity transmission spectrum even at finite thermal 
photon number. The system makes it possible to study further prominent collective phenomena such 
as superradiant decay of an inverted ensemble or the building of a narrowband stripline micromaser 
locked to an atomic hyperfine transition. 

PACS numbers: 37.30.+i, 42.50.Pq 



I. INTRODUCTION 

The idea of resonant coupling of an ensemble of atoms 
to a single cavity mode has been addressed in numer- 
ous aspects and contexts, some dating back several 
decades [l|. Recently, in the context of quantum infor- 
mation processing, such Hamiltonians attracted renewed 
attention because the ensemble can serve as quantum 
memory with long coherence times |2t5|- Despite small 
coupling of individual atoms, the strong collective cou- 
pling of the ensemble to a particular cavity mode allows 
for the coherent transfer of an excitation to the ensem- 
ble, its storage and its retrieval after some time shorter 
than the coherence time of the system. Hence, due to 
collective effects, one can utilize atomic transitions and 
geometries for which the strong coupling regime would 
not be accessible otherwise. 

As a particularly striking example, one can even envis- 
age the use of states that are very weakly coupled to the 
field, for example, an optically forbidden hyperfine tran- 
sition, which only couple to the field via magnetic dipole 
interaction. What makes this idea attractive and pos- 
sibly feasible with current technology is the fact that it 
should be possible to fabricate high-Q stripline waveguide 
cavities on the superconducting surface of a microchip, 
which confine the microwave mode to a very small effec- 
tive volume and to simultaneously trap a large ensemble 
of cold atoms very close to the surface. The combination 
of high-Q stripline waveguide cavities and atom-trapping 
technology surely will involve new challenges, but there 
seem to be no fundamental problems. As already demon- 
strated, such a cavity can be strongly coupled to on-chip 
Cooper-pair box qubits [6J. By combining the two sys- 
tems, one thus could establish a connection between the 
atomic ensemble and solid-state qubits. This setup hence 
bridges an enormous range of time scales starting from 
the sub-microsecond scale of solid-state qubits, over the 
millisecond lifetime of microwave photons, to the atomic 
hyperfine coherence lifetime of seconds. 



In the particular setup discussed here, the ensemble 
consists of a cloud of ultracold '^^Rb trapped in an on-chip 
magnetic wire trap and pumped to one of the trappable 
hyperfine levels, for example, F = \, mp = —1. The 
interaction between the atoms and the field is dominated 
by the magnetic dipole transitions between |_F = 1,771^) 
and \F = 2,m'p). These transitions are widely used for 
hyperfine manipulations of cold atomic ensembles by ex- 
ternally injected microwaves [7|. We assume in the fol- 
lowing that the experimental setup guarantees that the 
cavity is resonant with only one of the possible transitions 
(e.g. mp = — 1 and m'p ~ 1 with transition frequency 
uja/ (Stt) ~ 6.83 GHz, corresponding to T w 330 mK), 
and hence allows for the atoms to be treated as two-level 
systems. Actually, in some cases it is more favorable to 
use Raman-type coupling employing an extra radio- wave 
field to choose a suitable microwave transition Q . 

We ignore some of these technical details at this point 
and focus on the three main topics: After the introduc- 
tion of the model in Sec.|lll we first investigate conditions 
for strong coupling between the ensemble and the cavity 
and the experimental consequences when one adds the 
obscuring effects of thermal photons due to a finite cav- 
ity temperature. In Sees. IIII Al and lIHBI wc discuss the 
methods wc use. whereas in Sees. IIII Cj lIHFI wc address 
several aspects of the resulting dynamics. Here the op- 
tically aligned ensemble, which has much lower effective 
temperature, can be expected to act as a heat sink for 
the cavity mode removing thermal photons. As the up- 
per and lower hyperfine states have a virtually infinite 
lifetime compared with other system time scales, we can 
also completely invert the system, mimicking an effective 
negative temperature, and use it to pump energy into 
the system. As a prominent example, we study in Sec. lIVI 
the superradiant decay of a fully inverted ensemble again 
with some thermal photons initially present. Finally, we 
exhibit in Sec. |V] the possibility of building an ultra- 
narrow linewidth single-chip stripline micromaser oper- 
ating directly on an atomic clock type transition, which 
is in close analogy to an optical-lattice-based setup, as 



recently suggested in p. 



II. MODEL 

A. Collective atom-field Hamiltonian 

A single atom, formally represented here by a two- 
level system resonantly coupled to a cavity mode, can 
be well described by the Jaynes-Cummings Hamiltonian. 
For N two-level systems trapped so close to each other 
in the cavity that they see the same field and thus are 
coupled to the mode with equal strength g, wc then get 
the generalized Hamiltonian: 
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with a being the annihilation operator for a cavity pho- 
ton, (T^ being the excitation operator for the jth two- 
level system and [cj^if^] = '^l^ij- The frequency of the 
two-level systems and the mode are denoted by lOa and 
LOrm respectively. The coupling strength g ~ B (f) ■ fl/h 
depends on the strength of the magnetic field per pho- 
ton B at the position f of the atoms and the magnetic 
moment jl of the considered transition. 

What make an ensemble of atoms coupled to a cavity 
interesting are collective effects emerging from the com- 
mon coupling of all atoms to the same mode. This can 
be well illustrated by introducing collective atomic op- 
erators S — ^.jCrf and S^ = ^J^j'^j- The treatment 
in terms of collective operators provides a convenient ba- 
sis for classifying the possible states of the ensemble and 
is therefore discussed here. As wc see in Sec. IHBl we 
have to resort to Hamiltonian ([T]) in our particular treat- 
ment. The introduction of S and S^ leads to the Tavis- 
Cummings form of this Hamiltonian; 
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where single photons are coupled to distributed (delo- 
calized) excitations in the ensemble [l|. Let us shortly 
review some of its most known properties here. Math- 
ematically, the collective operators follow the standard 
commutation relations for angular momentum opera- 
tors S = (S^,S^S^), with S=^ = (S^±iS^). The cor- 
responding eigenstates of S^ and S^ are the so-called 
Dicke states |J,M), with S^\J,M) = J(J + 1)| J, M) 
and S^|J,M) = M\J,M), where J = 0,1,..., N/2 
and M = —J,...,J. Formally, a fully inverted ensem- 
ble corresponds to the maximum angular momentum of 
J = N/2 jl3,[n| and projection M = N/2. Repeated ap- 
plication of the collective downward ladder operator S^ 
on the initial state \J,J) = |e, e, . . . , e) gives the lowest 
state I J, -J) = \g,g,...,g). 

The states in between are generated according to 



S± I J, M) = ^( J ±M + 1){Jt M) I J, M±l) . (3) 



The interaction can then be conveniently rewritten in 



terms of normalized collective operators S = -7^ ^.^ a, 
to obtain 



± 



H = hajnia'a + hoJa S S -l-%off ( S a + a' S 



(4) 



with 5off = gyN. Note that in the case where the atoms 
in the ensemble couple to the cavity with different cou- 
pling constants gi, we generalize to S = j—.^igi<^i , 

with gctf = v^Ci5i ■ This reduces to geff = gvN if all gi 
are equal. To simplify matters, we remain with the case 
of equal coupling strength. 

Allowing only one excitation in the system, we see 
that the ground state |0)^ = \J,—J) = \g,g,...,g) 
is only coupled to the symmetric atomic exci- 
tation state S |0)^ = |1)^ = |J, -J+1) = 

-^{\e,9,9,---,9) + \9,e,g,...,g) ^ h \g, . . . ,g,e)), 

while other atomic states with only one excitation play 
no role. Hence, in this form we end up again with 
a two-level atomic system, where the dependence of 
the atom-cavity coupling on the number of atoms is 
explicitly visible. Even for transitions with a very small 
coupling constant g, strong coupling can be achieved for 
sufficiently large N. 

Note that in Eq. (g]) we use S « -i + ^ J , where 
S 



2 ' N 

2^ ^j erf . This approximation is exactly valid only 
either for a single atom or the special case where we con- 
sider only one excitation in the system. The constant — i 
is neglected in the Hamiltonian. In general, wc find for 
a state with J = y ^^'^ M — —J + s 



(J,-J + ,s|S \J,-J + s) = -- + — (5) 
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For s <S^ N wc neglect the second term on the right-hand 
side of Eq. ^ and find the approximation for S , which 
becomes exact for s = 1. 

For large ensembles with few excitations this approxi- 
mation is closely related to the bosonization procedure. 
For AI K, —J with J ~ N/2 and from 
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we find that for few excitations it is possible to iden- 
tify S and S with bosonic creation and annihilation 
operators. Hence, we end up with a system of coupled 
oscillators, for which a great deal of solution techniques 
exist. 

Let us now come back to the atom-field interaction 
[Eq. (HI] . It is well known that the eigenstates are coher- 
ent superpositions of the two previously introduced basis 



states, where the excitation is located either in the mode 
or in the ensemble. Let |0)^ and |1)^ = a^ |0)^ be the 
possible states of the mode and |0)^ and |l)jj be the en- 
semble states. With cua = iOm, the two eigenstates then 
read 

l+)=^(|l>a|0>,„ + |0)Jl)J , (8) 

|->=^(|l>a|0>„.-|0)Jl)J , (9) 

and as expected are separated by the energy difference 
2gee- Of course, the system possesses more states con- 
taining essentially one excitation quantum, but those are 
not directly coupled to the ground state if we consider 
only collective operators. The collective operators cou- 
ple states within one J manifold, like the previously 
discussed manifold with maximum angular momentum 
J = N/2 and M = —J ... J. Taking into account the 
manifolds of states with J < N/2, one can see that in 
general there is a large number of states describing an 
ensemble with n excitations. In the forthcoming calcula- 
tions including spontaneous emissions, such states with 
J < N/2 can be populated as well [12] ■ In addition, we 
also do not restrict the dynamics to a single excitation. 



B. Master equation including decoherence and 
thermal noise 



In any realistic implementation of the preceding model, 
coupling of the thermal environment to the field mode 
and the atoms is unavoidable. This generates several 
sources of noise and decoherence we have to address to 
be able to reliably describe the dynamics. Despite its 
high Q value, the microwave resonator still has a non- 
negligible finite lincwidth k = 1/r. In other words, a 
stored photon is likely to be lost from the cavity after 
the time r. Similarly, atomic excitations are assumed to 
decay with a rate that is, fortunately, in our case negligi- 
bly small in practice [13l |. However, we have to consider 
trap loss of atoms leaving the cavity mode, which gener- 
ates an effectively faster decay of the atomic excitation, 
denoted by the rate 7a. This can be to some extent con- 
trolled by a suitable choice of the trapping states and 
trap geometry. An additional and in general quite seri- 
ous source of noise are thermal photons that leak into 
the cavity. For an unperturbed cavity mode they lead to 
an average occupation number of 



n{uJm,T) = 
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where T denotes the temperature of the environment. In 
principle, such thermal photons are also present on the 
atomic transition and lead to a thermalization of the op- 
tically pumped atomic ensemble. Fortunately, the weak 
dipole moment of the atom renders this thermalization 
rate so slow that it can be ignored at the experimentally 



relevant time scales. In principle, even this rate could be 
collectively enhanced, but it largely addresses collective 
states only very weakly coupled to the cavity mode. 

Putting all these noise sources together, we can use 
standard quantum optical methods to derive a corre- 
sponding master equation for the reduced atom-cavity 
density matrix |l4|: 



^-i^[--] 



C[p] 



(11) 



with the Liouvillian 
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We assumed here that direct thermal excitations of the 
atoms can be neglected due to the weak coupling of the 
hyperfine transition to the environment. The only signif- 
icant infiux of thermal energy thus occurs via the cavity 
input-output couplers (mirrors). Note that the part of 
the Liouvillian describing spontaneous emission rcficcts 
the assumption that the atoms are coupled to N statis- 
tically independent reservoirs. The main reason for this 
treatment is that the decay rate ja summarizes the very 
small decay rate of atomic excitations and the loss rate of 
atoms from the trap. Since the loss of individual atoms 
from the trap is a noncollective process, the independent 
reservoirs assumption is advisable. This part of the Liou- 
villian cannot be written in terms of collective operators, 
and therefore it will not conserve J [la]. Therefore, states 
with J < N/2, including dark states, become accessible. 



III. SIGNATURES OF STRONG COUPLING 

A decisive first step toward applications of such sys- 
tem is the precise characterization and determination of 
their limits. In particular, the experimental confirma- 
tion of sufficiently strong atom-field coupling compared 
to the inherent decoherence processes is of vital impor- 
tance. This has to be seen in connection with extra lim- 
itations induced by thermal photons in the mode, which 
in contrast to optical setups play an important role in 
the microwave domain. We thus need reliable methods 
to determine the atom number, their effective coupling 
strength, and noise properties. In particular, we want to 
find the minimum temperature requirements that would 
make it possible to observe strong coupling. 



A. Numerical solution for small particle number 

To get some first qualitative understanding of finite T 
effects, we study the coupled atom-field dynamics in the 



regime of strong coupling under the influence of thermal 
photons based on the direct numerical solution of the 
master equation. Of course, here we have to resort to 
the limit of only a few atoms with increased coupling per 
particle. Nevertheless, at least the qualitative influence 
of thermal photons will become visible. For the practical 
implementation, we rely on the quantum optics toolbox 
for Matlab to explicitly calculate the dynamics of the 
density matrix [16|, which allows straightforward imple- 
mentation of the Hamiltonian in Eq. ^ formulated in 
terms of the collective operators. 

The cavity is pumped by a coherent microwave field 
with frequency cj; and strength 77, which in the frame ro- 
tating with uji is represented in the Hamiltonian by the 
additional term Up ~ ih [ija^ — ^* o) • From the station- 
ary solution, we then determine the steady-state pho- 
ton number in the cavity for different frequencies of the 
pump field to determine the central system resonances, 
where the pump frequency matches the eigenfrequencies 
LJrn i (7off of the couplcd ensemble-cavity system. At zero 
temperature and weak pumping, wc get the well-known 
vacuum Rabi splitting showing two distinct resonances 
separated by 2(7eff- With increasing temperature and 
number of thermal photons, these two peaks will get in- 
creasingly broadened and reside on a broad background. 
Figure [T{a) illustrates the effect. 

To compare the dynamics obtained from the restricted 
Tavis-Cummings Hamiltonian Htc in Eq. ^ with the 
dynamics of Hamiltonian H in Eq. ([1]), we compare the 
results for both cases in Fig. [TJb). Even for the rather 
small atom numbers chosen here, the difference and thus 
the influence of the nonsymmctric states is hardly visible 
in this observable. 

The peaks in the photon number in principle stay 
visible also for higher temperatures, but they start to 
broaden and finally vanish. Regardless, the detection 
on the thermal background gets technically more chal- 
lenging. The intracavity steady-state amplitude of the 
field shows a similar behavior and a determination of 
the splitting becomes increasingly impossible, despite the 
fact that phase-sensitive detection (homodyne) can help. 
An example is shown in Fig. [2] This effect is expected to 
be less important if the ensemble contains a large number 
of atoms. However, this regime is not accessible for direct 
numerical simulations and we have to develop alternative 
semianalytic approaches. 



Truncated cumulant expansion of collective 
observables 



To overcome the system size restrictions of a direct nu- 
merical solution, we now turn to an alternative approach 
that does not rely on the simulation of the dynamics of 
the whole density matrix. Instead, we derive a system of 
coupled differential equations for the expectation values 
of the relevant system variables. The inversion of atom i 
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FIG. 1. (a) Steady-state number of photons (a' a) in the 
pumped cavity for different detunings A™, = a;,„ — i^i and 
different temperatures, obtained from Hamiltonian Htc- The 
parameters chosen were k = 1, A'' = 2, <; = 3, 7a = 0.05, 77 = 
0.1. With increasing temperature the two peaks indicating 
strong coupling are superimposed by thermal photons. To 
compare the dynamics of Htc and H, we plot both results for 
T = 0.04 K (n = 3 • 10~*) in (b). 



then obeys 



Tr{a,f-p} 

-i2g ((a+a) - {cj- a})) j, {{at) + 1) (13) 



which couples to {<J^a) and (^g^ a^) ~ {'^t ^) ■ We as- 
sume that all atoms are equal, which allows us to replace 
(erf) with {al). Expectation values for pairs of different 
atoms like {o't'^j) '^^^ ^^ replaced with (cr^cr^). 

While these equations are exact in principle, the pro- 
cedure ultimately leads to an infinite set of coupled equa- 
tions. We thus have to start approximations and truncate 
this set at a chosen point, neglecting higher-order cumu- 
lants [^, |l3l ■ The truncation has to be carefully chosen 




FIG. 2. Effects of higher temperatures: For T = 0.1 K, (a) shows the steady-state number of photons, whereas (b) and (c) 
show the real and imaginary part of the steady state field in the cavity. The peaks in the photon number start to broaden 
and finally vanish. The amplitude of the field shows similar behavior. Panels (d), (e), and (f) show the same quantities for 
T = 0.7 K. The remaining parameters were chosen as in Fig. [T] 



and tested in general. Here we stop at third order, which 
in similar situations has proven to be well suited to de- 
scribe the essential correlations Q- 

The expansion for an expectation value of the form 
{ah) is the well-known relation {ah) = {ah) ^ -\- (a) (6), 
with {ab)^ being the covariance between a and h. Along 
this Une, one expands third-order terms in the form 



{ahc) = {abc)^^ + (a5)^ (c) + {ac) ^ {h) 

{a) {h) (c) . (14) 



neglected 



We make one exception in this expansion when it comes 
to the quantity i^a'^aal); the reason for this is discussed 
in Appendix [Bl 

The number of equations depends on the order of the 
cumulants we wish to keep track of. Furthermore, the 
problem is greatly simplified if there is no coherent input 
field driving our cavity. In this case, no defined phase 
exists in our system, so we can assume that (a) = (a^^) = 
(ct]^ ) = 0. Note that while for a single system trajec- 
tory a coherent field can build up as in a laser, for an 
average over many realizations the preceding assumption 
holds. For a covariance like {cr^a') = {cr^a') — (cr^) {a), 
we therefore find (cfi a\ = (cr^ a\ . The four remaining 



equations are 

d 
dt 
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+ {N-l){a+a^) 



^(ata) = -i.giV((aVr)-(aa+)) 



(16) 



—2k (a* a\ -\- 2Kn 



— (cr+o-2 ) = -la {crt(J2 ) 



(17) 



-ig(ai^)((arat)-(a+a)) . (18) 

To inject energy into our system without losing the prop- 
erty of having no defined phase, we can introduce an 
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FIG. 3. The steady-state field and photon number in the cavity for different frequencies of the pump. The size of the ensemble 
was chosen to be A'' = 10^ and we set k = 7 • 10^, 77 = 5 • 10^, 7a = 0.3 and g = 40. Panel (a) shows the steady-state number 
of photons, whereas (b) and (c) show the real and imaginary parts of the steady-state field in the cavity for T — 0.1 K, which 
corresponds to n = 0.04. The remaining figures show the same quantities for T = 0.7 K, corresponding to n = 1.67. The 
results show that for a sufficiently large number of atoms strong coupling remains observable despite of the presence of thermal 
photons. 



incoherent pump of the atoms. In essence, this gives 
an additional term in the Liouvillian very much resem- 
bhng spontaneous emission in opposite direction. For- 



mally, it reads -f EjLi \"j -j 
where w denotes the rate of the 
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in Eq. (fTS)) . 



Introducing a coherent pump leads to a larger set of 
13 coupled equations for the quantities (a), (af), (cj ), 
(aal)^, (aaf)^, (o-j^ct^)^, (a^a)^, (aa^)^, (a^a^)^, 

(erf (T^)^, (CTicrj")^, (crfcrl)^ and (a^aaf); for details, see 
Appendix [A] In this case we transform into a rotating 
frame with respect to the frequency of the pump coi . This 
results in A„j = ujm — ^i for the detuning of the cavity 
and Aa = Wa — w/ for the detuning of the atoms with 
respect to the pump frequency. 

In general, the set of equations is too complex for a di- 
rect analytical solution and has to be integrated numeri- 



cally. In this way we obtain the steady-state expectation 
values of relevant observables as the occupation number 
of the cavity (a^ a^ or the inversion of the ensemble {crl) ■ 
To compare the results of the obtained equations with 
the results in Sec. IIII A| wc plot the steady-state num- 
ber of photons and the field in the cavity for different 
frequencies of the pump laser in Fig. [31 We clearly see 
that effective strong coupling appears for a sufficiently 
large number of weakly coupled atoms and can stay vis- 
ible also at higher temperatures. Further discussion is 
given in Sec. lIIIEl In Fig. [5] we schematically depict the 
setup including the described loss and pump processes. 



C. Cavity output spectrum at finite temperature 

Naturally, the total field intensity in the cavity is only 
part of the story and significant physical information can 
be obtained from a spectral analysis of the transmitted 
field. Using the quantum regression theorem, the spec- 



trum of the light transmitted through one of the mirrors 
can be expressed in terms of the Fourier transform of the 
corresponding autocorrelation function of the field am- 
plitude. 



-^1W^ 





to the spectrum that arise from the reservoir, the cav- 
ity and cavity-reservoir interference. The initial condi- 
tions necessary for the Laplace transform are the steady- 
state values obtained either numerically for the coher- 
ently pumped cavity or analytically (see Sec. lV)) . 

We show the spectrum of the cavity without any pump, 
coherent or incoherent, but with T = 0.1 K in Fig. [51 The 
spectrum shows absorption dips at the frequencies of the 
coupled ensemble-cavity system. Some thermal photons 
that leak into the cavity are absorbed and lost into modes 
other than the cavity mode. In this form the thermal field 
is a broadband probe of resonant system absorption. 



D. Cooling the field mode ^vith the atomic 
ensemble 



FIG. 4. To simplify matters, we depict the cavity as a Fabry- 
Perot cavity which can be pumped through a mirror with 
high reflectivity. The observation of the dynamics is carried 
out using the second mirror, which has a lower reflectivity. 
Additionally, we can pump the ensemble incoherently from 
the side. 

At finite T this is not the full story, and to obtain 
the actual spectrum of the light impinging on the de- 
tector, one has to include the thermal photons in the 
output mode rcfiected from the cavity. Hence, the nor- 
malized first-order correlation function of the field out- 
side the cavity will have additional contributions from 
the correlation function of the thermal field, as well as 
of the correlation between the thermal field and the cav- 
ity field [18|. The latter is causing interference effects 
between cavity field and the thermal field. The corre- 
lations between reservoir operators and cavity operators 
can be expressed in terms of averages involving cavity 
operators alone [la, [l^. For a cavity radiating into a 
thermal reservoir, we find for the normalized first-order 
correlation function 
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with 
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(20) 



Here, r/ denotes the annihilation operator of a reser- 
voir photon and g {uj) denotes the density of states in 
the reservoir. Equations for the correlation functions 
in Eq. (|19p can be obtained via the quantum regres- 
sion theorem. The resulting system of coupled equa- 
tions is Laplace transformed to give the contributions 



The spectra depicted in Fig. [5] show a weak loss of 
thermal photons from the coupled ensemble-cavity sys- 
tem. Cavity photons are absorbed and sometimes scat- 
tered into a mode other than the cavity mode. As the 
ensemble can be nearly perfectly optically pumped into 
a particular state, its effective temperature is close to 
zero and hence well below the mode temperature. A 
relative purity of the ensemble of 10~^ corresponds to 
T - 28 mK, where we used huja/hT = In (lO^^) with 
uja/'^'^ = 6.83 GHz. This leads to the ciuestion as to 
what extent the thermal occupation of the mode can be 
reduced by thermal contact between the two systems via 
such energy transfer and loss. In Fig. [HJa) we show the 
dynamics of the photon number in the mode at different 
temperatures after putting the systems into contact. In 
Fig-EUb) we consider different loss or decay rates 7a of the 
excited atoms. In practice one could think of coupling to 
an magnetically untrapped atomic state or adding some 
repumping mechanism to increase this intrinsically very 
low rate. The dynamics is found numerically by inte- 
grating Eqs. (fT5 |l -p8 l) . To see the effect for increasing 
temperature in Fig. [5]Ja), we initialize the ensemble with 
all atoms in the ground state, whereas the mode con- 
tains n{uJm,T) photons. The decay rate of the atoms 
is chosen to be 7a = 5-10'*. With increasing tempera- 
ture the initial number of photons also increases. Due 
to coherent transfer and decay via the atoms, a constant 
fraction of the photons is removed from the cavity mode. 
In Fig. [6] (b) we show the same effect except that we now 
keep the temperature fixed to T = 4 K and vary the de- 
cay of the atoms 7a = 1 • 10'^ . . . 2 • 10^. The steady state 
of the photon number strongly depends on 7a. The red 
curve (with diamond markers) is the expected number of 
photons remaining in the cavity 



1 
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2k 



1 



(21) 



which coincides with the numerical results. The inversion 
(erf )gg can be calculated analytically from Eqs. ((T5)) - (fT5)) . 
The loss of thermal photons is proportional to the loss 
rate 7a and the number of atoms in the excited state 
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FIG. 5. (Color online) Overview of the transmitted spectrum S for different sizes of the ensemble (a). The temperature of 
the cavity is set to T = 0.1 K (n = 0.04). The remaining parameters were chosen to be k = 7-10'^, 7a = 0.3, g = 40, 
ijJa = i^m = 2-n ■ 6.83 GHz. The dips in the spectrum indicate that thermal photons are absorbed by the ensemble and re- 
emitted into modes other than the cavity mode. The increasing distance between the absorption dips reflects the increasing 
number of atoms. In panels (b)-(d) we depict the spectra at A^i = 3.2 ■ 10^, N2 = 1 ■ 10^, and N3 = 3.4 ■ 10^, indicated in (a) by 
the dashed horizontal lines. 
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The latter becomes very small if 7a be- 
comes large. Hence there is an optimal loss rate for each 
set of parameters. 

The removal of thermal photons becomes more effec- 
tive if the number of atoms is increased. Meanwhile, 
the inversion of the ensemble also drops since the frac- 
tion of excited atoms is decreased. Overall, the effect is 
clearly visible (sec Fig. [7]), but it seems that for the ac- 
tual parameters here its practical value remains limited. 
However, with a larger atom number and more tailored 
decay rates the method could be employed to reset a par- 
ticular mode shortly before starting any quantum gate 
operation. Note that this treatment of the cooling pro- 
cess is limited to short time scales since the permanent 
loss of excitations via 70 involves the loss of atoms from 
the ensemble. The number of lost atoms after the time 
i, approximated by the number of thermal photons that 
entered the cavity tnn, has to be much smaller than the 
ensemble size iV, which restricts the time t. 

In the situation where the mode is at T « OK and 
the atoms are subject to incoherent pumping we find in- 
creased transmission for the resonance frequencies (see 
Fig. [5]). We again recover the a/ZV dependence of the 
splitting of the peaks. 
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FIG. 7. Steady state of the photon number (a) and the in- 
version (b) for a loss rate 7a = k = 7-10'^. Increasing the 
number of atoms A'' leads to a more effective removal of ther- 
mal photons. The remaining parameters are chosen to be 
g = 40, uja=ujm = 2Tv 6.83 GHz, r = 4 K (n = 11.7). 



FIG. 6. (Color online) Loss of photons from the cavity mode. 
(a) Dynamics of the occupation of the mode for different 
temperatures and "/a ~ 5 ■ lO'' . A constant fraction of the 
photons is removed from the cavity mode, (b) For fixed 
T — 4K {n = 11.7), the loss rate of the atoms is varied 
between 7a = 1 • 10^ and 2 ■ 10^ . The steady-state number 
of photons shows that there is an optimal loss rate. The 
red curve (with diamond markers) corresponds to the ex- 
pected number of photons according to Eq. (|21|l . The re- 
maining parameters were chosen to be k = 7 • 10^, g = 40, 
Wo = a 



27r- 6.83 GHz, N = 10^ 



10 




0.5 0.5 1 X 10^ 

Lc(rad/s] 



FIG. 8. (Color online) The temperature of the cavity is now set to T = 0.001 K (n = 0) and an incoherent pump of the atoms 
with w = 0.05 is switched on. The spectrum now shows increased emission at the frequencies of the coupled system. Again 
figures (b)-(d) depict the spectra at A'^i = 3.2 • 10^, A'^2 = 1 • lO'', and A'^3 = 3.4 • 10^, indicated in (a) by the dashed horizontal 
lines. 
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FIG. 9. Steady state field in the driven cavity, real part (a) 
and imaginary part (b). The size of the ensemble is chosen to 
be A^ = 10^ The fines for T = 0.01 K (n = 0) and T = 10 K 
(n = 30) coincide. 



E. Coherently driven cavity mode 



FIG. 10. Steady state field in the driven cavity, real part 
(a) and imaginary part (b), for T = 0.01 K (n = 0) (solid 
lines) and T = 10 K (n = 30) (dashed lines). For the small 
ensemble with A'^ = 10^ atoms we recover the effects of the 
thermal photons. To compensate for the lower number of 
atoms, the coupling is chosen to g = 1200. Otherwise, the 
splitting would be covered by the cavity linewidth. 



An experimentally readily accessible quantity is the 
cavity field amplitude, which can be deduced by phase- 
sensitive (homodync) detection of the output. This quan- 
tity is much less obscured by random thermal field fluctu- 
ations than the spectral intensity in total. As a phase ref- 
erence, we therefore now introduce a coherent phase sta- 
ble pump of the cavity, which is again represented in the 
Hamiltonian by the additional term Hp = ih (^rja^ — rj*a). 

As mentioned previously, a coherent pump strongly in- 
creases the number of nonvanishing cumulants and at our 
level of truncation leads to a set of 13 coupled equations, 
which can be found in Appendix |Al Based on this set, 
we can calculate the stationary real- and imaginary part 
of the field in the cavity after transient dynamics. The 
amplitude of the field inside the cavity becomes maximal 
if the frequency uji of the driving laser hits one of the 
resonances of the coupled system. As we give a phase 
reference now the effect of a higher temperature on the 
field in the cavity is barely visible, in particular if we 
chose a large ensemble of A^ = 10'' atoms (see Fig. [9]). 

Note that although not giving the vacuum Rabi split- 
ting, the average atom-field coupling can be still deduced 



from these resonances as g enters in their frequency. If 
we go back to a rather small ensemble of A^ = 10^ atoms, 
the influence of the temperature becomes visible. To en- 
sure that we still can observe well split levels, which are 
not covered by the linewidth of the cavity, we increase 
the coupling constant g in our simulation. The results 
in Fig. [To] show that thermal effects become visible in 
the field only if the number of thermal photons is not 
negligible compared to A^. 



F. Spectrum of the coherently driven cavity 

The spectral intensity distribution of the coherently 
pumped cavity is calculated in a way similar to that de- 
scribed in Sec. IIII CI In contrast to the incoherent pump 
process that excites atoms in a noncoUective way, as can 
be seen from the Liouvillian, the interaction with the 
coherently pumped mode is a collective interaction. To 
demonstrate this behavior, we calculate the spectral dis- 
tribution of the mode intensity, without caring about the 
reservoir it radiates into, and the spectrum of the fluo- 
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FIG. 11. (Color online) Incoherent spectrum of the mode (a) and of the fluorescence of the ensemble (b) for T = 0.025 K, 
K = 7 ■ W^, 7 = 0.3, N = 10^, 1] = 9 ■ 10^ in red (gray) and r; = lO** in black. The pump driving the system is on resonance 
with the cavity and the ensemble. Insets (c) and (d) show a magnification of the central peak of each spectrum. The lower red 
(gray) line represents the result for ?) = 9 ■ 10^ . 



rcsccncc of the atoms. 

The cavity mode and the atomic transition are assumed 
to be on resonance with uja = ^m = 27r • 6.83 GHz, which 
is also the frequency of the pump laser. To calculate 
the incoherent part of the spectra we need the Fourier 
transform of the two-time correlation functions 



lim (a) {t)a{t + T)) = 

lim {{a^ (t) ait + r)) - (a^t)) {a [t + r))) (22) 



and 



^lim(a+(i)a-(t + r))^ = 
lim ((a+ (t) a- {t + r)) - (a+(t)) {aj {t + r))) 



(23) 



The quantum regression theorem and Eqs. (jAip to (jASp 
give 



-(at(0)a(r))^ = 



At 



{k + iA„) (at (0) a (r))^ - \gN {a^ (0) aT (r))^ (24) 



and 



dr 



(c7+(0)a-(r))^ 



la 



-(f +iA,j(a+(0)a7(r))^ 
ig ({a{T)) {a+ (0) a] (r))^ + (^.f (r)) {a+ (O)a(r)), 



(25) 



lim. 



{aHt)ait + T))^ 



couple to four 



where we us 
(at (0) a (r))^ and Hmt^oo {(^t (*) <^j 
(a+ (0) ay (r)) . Equations ^ and dH 
other two-time correlation functions that have to be cal- 
culated. To solve for the desired quantities, we Laplace 
transform both sets of equations and use Cramer's rule 

to obtain (at (0) a (t))^ (s) and {a+ (0) aj (r))^ (s). The 
necessary steady-state values are obtained numerically. 

The incoherent spectra of the mode and of the atoms 
both show a narrow peak at uja = ujm which has a width 
of sa 27a (sec Fig. [TT|) . The double-peaked structure is a 
remainder of thermal excitations acting as a broad band 
probe for the ensemble-cavity system. With increasing 
strength of the coherent pump the narrow central peak 
becomes dominant. The appearance of the central peak 
is probably related to weak contributions from almost- 
dark states (very weakly coupled to the mode). Let us 
mention in this context that the atomic ensemble is not 
restricted to a manifold of the Dickc states with fixed J 
since we include spontaneous emission in our model. It 
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FIG. 12. (Color online) Dynamics of superradiant emission: 
numerical solutions of the dynamical equations for an ensem- 
ble of A*' = 10^ atoms. The rapid drop of the inversion {ai) 
during the emission can be seen in (a). The exchange of exci- 
tations between the ensemble and the cavity is characterized 
by {(Kyf) in (b). Negative imaginary part of {ao-t} indicates 
emission from the ensemble into the cavity, where a positive 
imaginary part indicates absorption of cavity photons by the 
ensemble. In (c) the number of photons in the cavity is de- 
picted. Panel (d) shows the spin-spin correlation (o"^o"2^). In 
this example the temperature of the mode was chosen to be 
r = 4K. 



is hence possible that the ensemble ends up in a dark 
state, where it does not couple to the cavity mode. The 
time-constant that determines the decay into and the de- 
cay of such a dark state is of the order l/ja- This allows 
for the buildup of long time coherences, and the times 
the ensemble is in a dark state significantly change the 
statistics of the photon emission. The result is then a 
narrow peak in the incoherent spectrum J19l |20[ , where 
the width of the peak is determined by the characteristic 
time the ensemble resides in a bright or dark state, in 
our case ja- 

In the case of an incoherently pumped ensemble, the nar- 
row peak does not arise. A reason for this can be the na- 
ture of the incoherent pump which is noncollective and 
hence able to pump the ensemble out of a dark state in 
a shorter time. This is not possible in the case of the 
coherently pumped cavity: Spontaneous emission brings 
the ensemble to a dark state, but the collective interac- 
tion with the mode cannot reach it there. 
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FIG. 13. Dynamics of the photon number in the cavity with 
increasing temperature. The onset of the superradiant emis- 
sion is shifted to earlier times since the initially present ther- 
mal photons contribute to the fluctuations that trigger the 
emission process. 



pare an almost fully inverted ensemble, which can feed 
energy into the system and corresponds to an effective 
negative temperature [2l| . Since we have no initial phase 
bias in the system, Eqs. p5|) -p ^ are suitable for de- 
scribing the dynamics. The resulting superradiant dy- 
namics for an initially fully inverted ensemble is depicted 
in Fig. [121 In free space the emission occurs in a charac- 
teristic burst of duration « —^ [221 . The presence of the 
cavity causes a partial reabsorption of the emitted pho- 
tons, which can be seen in Fig. [T^ (c). Due to the large 
number of emitted photons it should be clearly detectable 
even on a fairly high thermal background. Following the 
pulse shape, one also can extract the effective coupling 
parameters to characterize the system. 

The process of superradiance can create a transient 
entangled state of the ensemble [23[ . This entanglement 
can be revealed by entanglement witnesses which can be 
inferred from the calculated observables. We have seen 
some indication of such entanglement appearing. How- 
ever, the persistence of the entanglement under the in- 
fluence of noise and with the presence of the cavity will 
be part of future work. 

The onset of superradiant emission is determined by 
spontaneously emitted photons that trigger the forth- 
coming burst of radiation. The presence of thermal pho- 
tons is expected to reduce the time until the onset of the 
burst. This behavior is recovered by our equations as 
shown in Fig. 1131 where we depict the dynamics of the 
photon number in the cavity for different temperatures. 



IV. SUPERRADIANCE 

A great advantage of the considered setup is that one 
has full control of the atomic state. Hence, instead of 
starting at a zero-temperature ground state we can prc- 



V. NARROWBANDWIDTH HYPERFINE 
MICROMASER 

The collectively coupled ensemble can be used to con- 
struct a stripline micromaser with a very low linewidth. 
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FIG. 14. (Color online) Steady-state inversion of the ensemble (a) and occupation of the cavity mode (b) for varying ensemble 
size A^ and pump strength w. The temperature was chosen to be Ti — 0.001 K, which corresponds to an empty cavity {n = 0). 
Vertical red dashed lines mark the masing threshold w = 'ya ~ 0.3. The horizontal red dashed lines at A'^ = 10^ indicate the 
position of the curves shown in (c) and (d). From (c) we recover the passage of the inversion through zero for w — 'ya ~ 0.3. 
At this point we see in (d) a rapid increase of the photon number in the mode. 



To this aim the inversion of the ensemble is sustained by 
an external incoherent pump of the atoms. In contrast to 
the calculations in Sec. IIIIC) wc ignore the fact that the 
cavity radiates into a thermally occupied reservoir. Af- 
ter passing the masing threshold, the thermal occupation 
outside becomes negligible. To determine the linewidth 
of the emitted light we calculate the Laplace transform 
of the two-time correlation function (a'''(t)a(0)). Using 
the quantum regression theorem we find 






-1.9^^1/ 



igN 

ia-\--ya. 



(at(i)a(O)) 
(a+(t)a(0)) 



(26) 



Laplace Transform of Eq. (j26p gives 



-igN 



15 (^i; 



W+Ja 



+ s 



,(a+(iMO)), 






(27) 



where {■)gg denotes steady state values and ~ denotes 
Laplace transformed quantities. 

The steady-state values on the right-hand side of 
Eq. (|27p can be obtained analytically. Setting the 
time derivatives of the dynamical equations to zero, a 



quadratic equation for (cf )^^ is attained. One of the solu- 
tions yields a physically meaningful result for calculating 
the remaining steady-state values and hence (a^a) and 

{(^lo) ■ To illustrate the effect of the increasing pump 
strength, we show the steady-state inversion of the en- 
semble and the occupation of the cavity in Figs. [14] and 
1171 Once the critical pump strength is reached, the sys- 
tems behave identically for different temperatures. The 
number of atoms in the ensemble is varied between 10'^ 
and 10^, whereas the pump parameter w ranges from 
10"'^ to 10^. In both figures we mark the pump strength 
w — "/a ~ 0.3, for which we find the inversion becomes 
positive, with a vertical red dashed line. At this point 
we also find a rapid increase of the number of photons in 
the cavity. The horizontal lines mark the cross sections 
for N = 10^ shown in^what follows. 

To solve for (a'l'(i)a(O)), we use Cramers rule on 
Eq. dSll), which yields 



(at(i)a(0))(.) = 






(28) 



so that with s — —iu the spectrum is given by 



S(w) = 



1 

2^ 



;at(i)a(0))(..) + (at(i)a(0)) (..) . (29) 
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FIG. 15. (Color online) (a) Linewidth of the spectrum S(a;). 
For each set of parameters we numerically determine the 
linewidth of the spectrum. The parameters were chosen to 
be K = 7■10^ 7a = 0.3, g = 40, Wm ^ c^a = 27^•6.83•10^ 
T = 0.001. (b) Exemplary spectrum for N = lO'' and 
w — 0.55, marked in (a) by the white cross. 



FIG. 16. (Color online) Finite temperature effects in the spec- 
trum. Panel (a) again shows the linewidth of the spectrum. 
Below the critical pump strength we recover the linewidth of 
the cavity. The parameters were chosen to be k = 7 ■ 10^ , 
7a = 0.3, g = 40, ujg = uja ^ 27r • 6.83 • lO'', T = 0.1. (b): 
Exemplary spectrum for A*' = lO'' and w — 0.55, marked in 
(a) by the white cross. 



For each set of parameters we calculate the spectrum 
and determine the linewidth numerically. The linewidth 
of the maser for two different temperatures Ti = 0.001 K 
and T2 = 0.1 K is shown in Figs. [TSl a) and [TBT a). For 
w = 7a = 0.3 we see a rapid drop in the linewidth for 
both temperatures and a resulting minimal linewidth of 
S ~ ^4.7 • 10"'^ Hz. Above the critical pump strength 
the pump noise destroys the coherence between the in- 
dividual atoms 0. In in Figs. [TSTb) and [EUb) we plot 
exemplary spectra for N ~ 10^ and w — 0.55, marked by 
the white cross. 

For T2 = 0.1 K the cavity contains on average h = 
0.04 photons that can be recovered from the constant 
background in Figs. ITTT b) and [TTTd). Since the number 
of thermal photons is small compared to the considered 
ensembles, the inversion is nonsensitive to the increased 



temperature. 

In Fig. [16] we recover the linewidth of the cavity 
K = 7 • 10^ if the pump is below threshold and again if 
the pump exceeds a critical strength Wmax- Above Wmax 
the coherence between different spins is destroyed by the 
pump noise [9| . The behavior between the threshold and 
Wmax resembles the behavior for T — 0.001 K shown in 
Fig.[l5l 



VI. CONCLUSIONS 

Our studies show that a hybrid cavity QED system 
consisting of a stripline microwave resonator at finite T 
and an ensemble of ultracold atoms is a rich and versa- 
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FIG. 17. (Color online) Same as Fig. [T4]witli T2 — 0.1, which corresponds to n = 0.04. In (d) the thermal excitations in the 
mode appear as a constant background from which the increase due to the pump stands out. 



tile setup for observing and testing prominent quantum 
physics phenomena. The effectively very cold tempera- 
ture and good localization of the atomic cloud allow sym- 
metric collective strong coupling to the microwave mode. 
While the weak magnetic dipole coupling requires large 
atom numbers and an extremely well localized microwave 
mode to obtain significant coupling, it also makes the sys- 
tem quite immune to external noise. In addition to the 
long lifetime of the atomic states, this renders the sys- 
tem an ideal quantum memory or allows for very narrow 
spectral response or gain. As all the atoms are identical 
and well trapped, the system exhibits only a very narrow 
inhomogeneous broadening. Operated in an active way, 
one thus can envisage a truly microscopic maser with an 
very narrow linewidth directly locked to an atomic clock 
transition. The uniform coupling and the possibility of 
efficient optical pumping enables the study of superradi- 
ant decay into the stripline mode, where a precise phase 
and intensity analysis of the emitted radiation can be 
performed. 

While many of our considerations are guided by param- 
eters expected from an ultracold atom ensemble, it is easy 
to generalize to alternative setups using NV-centers or 
other solid-state ensembles. There larger ensembles can 
be easily envisaged but one also gets much more varying 
coupling constants and larger inhomogeneous widths. It 
is not obvious whether the technically more simple setup 
and larger numbers in this case can compensate for these 
imperfections. 



This could be particularly important for a next step: 
possible optical readout of the ensemble. For the atomic 
case, the uniformity of the coupling over many optical 
wavelengths should allow a nice directional readout of 
the ensemble state, once a laser could coupled in. 
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Appendix A: Coupled Equations detuning of the cavity and Aq = Ua — uji for the detuning 

of the atoms. The coupled equations are given by: 

Transformation to a rotating frame with respect to the 
frequency of the pump ui results in A,„ = a;„i — ^i for the 

I 

— (a) = - (k + iA„0 (a) - igN (a^) + ry (Al) 

^ (aO =-i^- iA™) (aO + i.9^ i^t) + V* (A2) 
A ^al) = -2i, (((a+a)^ + (.+ ) (a)) - ((.faOc + (^D («0)) 

-7a(l + K)) (A3) 

^ {^i) = - ( Y + iAa) (^r) + ■^9{K<^)c + K) («)) (A4) 



dt 



(-1+) = - (y - iAa) (^n - i5 ((<^f at)^ + K) (at)) (A5) 



A ((aa+)^ + (a) (a+)) ^ _ (« + :| + i (A,„ - A,)) ((aa+)^ + (a) (a+)) 

- i.9 (^^i^ + (^ - 1) {{<rt-^)^ + {<rt) (-T.-))^ 
-ig(ajata)+77(a+) (A6) 



^ ((aaO, + K) (a)) = - (a. + iA^) ((aaj), + K) (a)) + r; K) 

- 27a ((a(T+)^ (ar) + (aa^)^, {<J+) + (a) {a+a^)) 

- ig [2 (2 {aU), («) + (aa)c K) + («) («> (^0) 

-(l + 2(ata)J(ar)-2(aar),(at) 

-2(ata-)^(a)-2(at)(a)(ar) 
+ (iV-l)((aJa-)^+K)(a-))] 



(A7) 



-7a(K'T2-)^ + (a+)K)) (A8) 



A ((ata)^ + (at) („)) = _ i,^. (((„t,-)^ + („t) (,-)) _ ((«,+ )^ + (,) (,+ ))) 

- 2k ((fltfl)^ + (flt) (a)) + 2k7i + 77* (a) + 77 (a^) (A9) 



^ ((«'^r), + («) ('^r)) = - (i (Am + AJ + K + |i) ((aar)^ + (a) (^rf)) 
- ig(iV - 1) ((^1-^3-)^ + (af) (af)) + 77 (ar) 
+ ig (2 {ala)^ {a) + {aa)^ (aj) + (a) (a) {al)) (AlO) 
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^ ((atat)^ - (at) (at)) = - {2k - i2A„0 ((atat)^ + (at) (at)) 

+ 2igN ((a+at)^ + (at) (a+)) + 2rf (at) 



(All) 



+ i2.g {{ala^)^ (a) + (aja)^ (af) + {a^ a) ^ K) + K) {'^l) («)) 



(A12) 



di 



((^i4). + K> (-0) = (- Y + iA.) ((aiV+), + K) (a+)) 



27a ((^^(Tf) (a+) + {ata+)^ (af) + (afal)^ (a+)) 

iff [2 ((^^2+), («> + 2 (^i+a)^ (a+) + (a+) (a+) (a) - (ara+)^ (at) 

K«0. (-n - K«0. (-D - K) (-D («0) + 

{alal)^^ (at) + 2 (ajat)^ K> + K> K> («^)] 



(A13) 



di 



(K^2)c + K) (^0) = -4i.9 (('^^^2+)^ (a) + {ata)^ (a+) + (a+a)^ K) + (a) K> ('^0 

- (a^a^,)^ (at) - (ajat)^ (a^) (a^a^)^ K> - («^) {^i) K)) 

- 47a ((^T+crr) K) + (a+a|)^ (af) + (a^crl)^ (a+)) 



(A14) 



Appendix B: Validity of the cumulant expansion 

The validity of the truncation of the expansion per- 
formed previously relies on the assumption that higher- 
order cumulants are negligible. This can be checked in 
principle by truncating at higher orders and comparing 
the results. In general, it turns out that there is one 
cumulant which requires more care: the correlation be- 
tween the inversion and the photon number (atacrj). In 



I 

the regime where ^^ > 1 and ^2att > i holds, the num- 
ber of photons necessary to saturate the ensemble is low. 
Therefore small fluctuations of the photon number can 
cause significant changes in the inversion |24l |. The cor- 
relation between the photon number and the inversion 
(^a''^a(jl'^ is therefore kept in our calculations. Hence, an 
expansion like in Eq. ([T4|) would not be advantageous be- 
cause none of the terms could be dropped. We therefore 
calculate the dynamical equation for (atacrj^) which gives 



— (atatJi) = (-2k - 7a) (ataffi) - 7a (flta) + 2Kn (al) 

— \g {^(^aai^ — (atcr^) + 2 (i^a^ aaaf^ — (atatacrf ))) 
-ig{N- 1) {{a^ala^) - {a<jla+)) 



(Bl) 



The expectation values of products of three operators 
are expanded as in Eq. ([T^ . except for (atacrf). The 
expansion of expectation values with four operators is 
more involved and produces also expectation values of 
products of three operators which are again expanded. 



I 

Cumulants of order three and four are neglected. The 
resulting equation for (atacrf ) can be integrated numer- 
ically along with the equations for the other quantities 
mentioned previously. 
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4 X 10^ 




-wilhuut (aUi(T-\ 
— with lu^uo-') 



2.5 3 X 10-^ 



To estimate the influence of the correlation between 
the photon number and the inversion i^a'aaf) on the 
dynamics we plot the photon number in the cavity dur- 
ing the decay of a fully inverted ensemble. We therefore 
integrate a set of 12 equations that is obtained if (a^acrf ) 
is expanded and (a^aaf') is neglected. For comparison 
we also show the dynamics obtained from the full set of 
13 equations in which ^a^atrj) is kept [sec Fig. [T51 (a)]. 



2 X 10 




FIG. 18. (a) Dynamics of the photon number in the cavity 
mode for an initially fully inverted ensemble. The solid line 
shows the dynamics produced by the set of equations where 
(a^ aal) was neglected. The dashed curve is the result of the 
full set of 13 equations in which (a^acrf ) is kept, (b) Numer- 
ically obtained cumulant (^a^aal^ in the steady state. With 
increasing loss rate of the cavity k the correlation between 
photon number and inversion decreases. 



The steady state of both solutions differs only slightly. 
The correlation (^a^aul) is shown in Fig. [18] (b). With 
increasing cavity decay rate k the correlation between 
photon number and inversion decreases. 
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